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Abstract 

We investigate the non-equilibrium charge and spin-dependent currents in a quantum ring with 
a Rashba spin orbit interaction (SOI) driven by two asymmetric picosecond electromagnetic pulses. 
The equilibrium persistent charge and persistent spin-dependent currents are investigated as well. 
It is shown that the dynamical charge and the dynamical spin-dependent currents vary smoothly 
with a static external magnetic flux and the SOI provides a SU(2) effective flux that changes 
the phases of the dynamic charge and the dynamic spin-dependent currents. The period of the 
oscillation of the total charge current with the delay time between the pulses is larger in a quantum 
ring with a larger radius. The parameters of the pulse fields control to a certain extent the 
total charge and the total spin-dependent currents. The calculations are applicable to nano-meter 
rings fabricated in heterojuctions of III-V and II- VI semiconductors containing several hundreds 
electrons. 

PACS numbers: 78.67.-n, 71.70.Ej, 42.65.Re, 72.25.Fe 



I. INTRODUCTION 



Study of the spin-orbit interaction (SOI) in semiconductor low dimensional structures 
and its application for spintronics devices have attracted much attention recently There 
are two important kinds of SOI in conventional semiconductors: one is the Dresselhaus 

n 

SOI induced by bulk inversion asymmetry |2|, and the other is the Rashba SOI caused by 
structure inversion asymmetry [3|. As pointed out in [4|, the Rashba SOI is dominant in a 
narrow gap semiconductor system and the strength of the Rashba SOI can be tuned by an 



external gate voltage in HgTe [5fl, InAs 



6j, In x .Gai_ x As [7], and GaAs 3, 8] quantum wells. 



Recent research is focused on the electrically-induced generation of a spin- dependent current 
(SC) mediated by SOI-type mechanism, e.g. as in the intrinsic s pin Hall effect in a 3D p- 
doped semiconductor [9| and in a 2D electron gas with Rashba SOI 10j. Here we study a high 



quality spin-interacting quantum rings (QRs) with a radius on the nanometer scale [ill. Il2|. 
These systems show Aharonov-Bohm-type (AB) spin-interferences [13I . 11^ . In particular 
we investigate the dynamics triggered by time- dependent electric fields as provided by time- 
asymmetric pulses [l5] or tailored laser pulse sequences jig . The quantity under study is the 



spin-resolved pulse-driven current, in analogy to the spin-independent case [17|, [18|, [19|] . In a 
previous work 20j, we investigated the dynamical response of the charge polarization to the 
pulse application. No net charge or spin-dependent current is generated because the clock- 
wise and anti-clock-wise symmetry of the carrier is not broken by one pulse or a series of 
pulses having the same linear polarization axis. This symmetry is lifted if two time-delayed 
pulses with non-collinear polarization axes are applied 17|. However, to our knowledge all 
previous studies on light-induced currents in quantum rings did not consider the coupling of 
the spin to the orbital motion (and hence to the light field), which is addressed in this work. 
As detailed below, having done that, it is possible to control dynamically the spin- dependent 
current in a ID quantum ring with Rashba SOI by using two time-delayed linearly polarized 
electromagnetic pulses. For transparent interpretation of the results only the Rashba SOI 
is considered in this work. The presence of the Dresselhaus SOI may change qualitatively 
the results presented here for the spin-dependent non-equilibrium dynamic of the carriers, 
which can be anticipated from the findings on for the equilibrium case 21 1. 
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II. THEORETICAL MODEL 



We study the response of charges and spins confined in a one- dimensional (ID) ballistic 



QR with SOI to the ar 
electromagnetic pulses [17 
Hi(t) Q, with 

,2 



ication of two short time-delayed linearly polarized asymmetric 
22]. The effective single particle Hamiltonian reads H' = -^goi + 



#soi = ^ + ^( r ) + x (a x p)z ' &lit) = ~ er ' E(t) + /iBB(t) ' & - 



E(£) and B(£) are the electric and the magnetic fields o: 



24 



;he pulse. Integrating out the r 



25 



26 



27] 



dependence -ffgoi rea ds in cylindrical coordinates |23|, 

2 0o 2dJ 2cj Ct'o 

9o, = 4>o = h/eis the flux unit, = i?7ra 2 is the magnetic flux threading the ring, a is the 



H, 



SOI 



H cr 2 



(2) 



radius of the ring, friu = h 2 /{m*a 2 ) = 2E Q , hu R = 2an/a, fajj B = 2\i B B and B are due to a 
possible external static magnetic field B = Be z . The single-particle eigenstates of -f^goi are 
represented as \£f(</?) = e^ n+1//2 ^z/ 5 (7, y>) where v s {^,ip) = (a 5 e~ l</,//2 , b s e ltp ^ 2 ) T are spinors 
in the angle dependent local frame, and = cos(7/2), h 1 ' = sin(7/2), = — sin(7/2), = 
cos(7/2), (T means transposed) where tan7 = —Qr = —^r/^q (if we ignore the Zeeman 



splitting caused by the static magnetic field |27J, l28|). 7 describes the direction of the spin 



quantization axis^ as illustrated in Fig. (la). The energy spectrum of the QR with the SOI 
reads Q, 



24 



25 



26 



23, 



Ei 



TlJjJn 



[n — 



+ 



1-Sw, 



Q 2 



R 



(3) 



w 



1 + Qr = 1/ COS 7, 

where S = +1 (S = —1) stands for spin up (spin down) in the local frame. 



III. PULSE-DRIVEN SINGLE-PARTICLE DYNAMICS 

We apply two time-asymmetric pulses to the system (see Fig. (lb)). The first one (at 
t — 0) propagates in the z direction and has a duration 7^. Its E-field is along the x direction. 
Td is chosen much shorter than the ballistic time of the carriers in which case the QR states 
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FIG. 1: (Color online) (a) Schematic graph of the geometry, spin configuration and the applied 
pulses is shown, (b) Time-delayed asymmetric pulses are schematically drawn, (c) Energy spectrum 
for a ring with spin orbit interaction. As defines the distance between the spectrum symmetry 
axis and the smallest nearest integer. 



develop as [17|, |22|, |29] 

^(ip, t>0) = $S(<p, t < 0)e iai cosv , ai = eap/h } p= - [* E(t)dt, (4) 

Jo 

where E(t) = Ff(t), F and f(t) describe the amplitude and the time dependence of the 
electric field of the pulse respectively. In the following, we use F% and F2 to characterize the 
first and the second pulses. The pulse effect is encapsulated entirely in the action parameter 
a%. With the initial conditions n(t < 0) = n and S(t < 0) = So and using Eq. (@J one finds 

*5(<M) ^E^K^o.tle^^e-'^/V 5 ), (5) 



(IS 



with 

C S = \ S SS $nn for t < 0, 

[ W n °" n ^o-n(ai) fort>0, 
where J n is the n-th order Bessel function. For the time-dependent energy we find 

B&(t>0) = E*{t<0) + ^, (7) 



with E%°(t < 0) is given by Eq. Applying a second pulse at t = r with the same duration 
Td but the electric field being along the y axis (see Fig. (lb)), the wave functions develop 
as i&n°((p, t > t) = ^^(y, t < r )e iCV2Smv: ', where a 2 is the action parameter associated with 
the second pulse. ^f^°((p,t = r~) follows from Eq. (jSj). For t > r the expansion coefficients 
behave as C n s /(n , S , t > r) = E n W^^no-nM Jn'-n(a 2 )}e^ E n'- E ^ h . 



IV. NONEQUILIBRIUM SPIN AND CHARGE CURRENTS 



A single pulse does not generate in QR any net charge current because of the degeneracy 
of the orbital states. However, the charge will be polarized 



20 



22l | and corresponding dipole 



moments oscillate in the x direction with an associated optical emission. Applying a second 
pulse as described above leads to a non-equilibrium net current, in addition to the persistent 
charge current caused by the static flux and the SOI which causes a SU(2) vector potential 



and manifests itself in an induced spin-dependent persistent charge current 
Consequently, a non-equilibrium spin- dependent current is induced. 



21, 



28 
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The line velocity operator is 



32] 



%h _ h 



. OCR 



m*a m*a 4>q h 



which is associated to the operator of the angular ve 



ocity v^/a [33]. Contributions to the 



persistent charge current from each QR level read 34j 

i r 27r 

In ,S 



*Pl drf no<So (r',t>r), 



ri 



where 

jS, l flb = ^[*5 ,t (^*)^*5( r '.*)]- 

Upon algebraic manipulations we find 

K,s (t > t) = I { Z So (t >r) + lV So (t > r) 



(8) 



The index "(0)" stands for the static persistent charge current (PCC) which exists in the 
absence of pulse field, whereas the index "(1)" indicates the pulse-induced dynamic charge 



current (DCC) 



351 ] without 



34 



'he 

m 



CC is caused by a magnetic U(l) flux and has been studied extensively 



37( ] or with the spin interactions 33]. It has been experimentally 



observed both in gold rings of radius with 1.2 and 2.0 /zm 



381 ] and in a GaAs-AlGaAs ring 



of radius about 1 /xm [39]. The SOI scattering effects were also studied [40]. The PCC 
carried by the states characterized by uq and So reads (please note the current in this work 
is defined as flow of positive charges, which is opposite to the direction of flow of electrons) 

r(o) \ — a r„L_l, 1 ~ S ° w 



C*(* > r ) = M° ( n ° - X + — ) > ( Q ) 



'0 



where J = 2E a/(j) is the unit of CC, the second term on the right hand side of Eq. (J9]) 
stems from the static magnetic field; the third term is a consequence of the SU(2) flux of 



the SOI [22! . The DCC P art is 

Ino,5o(* > r) = e^/o {a 2 (cos^)5(r)} , (10) 



where 



(J) X *S* TV 

(cos</?)f°(r) = «i/i(fi T i) sin6 r cos[2(n - — H )fe T ], 

b T = uj t/2, 



Or! = «iV2(l -cos(26 T )), 

/l(fi T l) = Jo(firl) + ^(firl). 

To obtain the total persistent charge current and the dynamic current we have to consider 
the spin-resolved occupations of the single particle states. For simplicity we operate at zero 
temperatures and ignore the relaxation caused by phonons or other mechanisms, i.e. we 
confine ourself to times shorter than the relaxation time. The general case can be developed 



along the line of Ref. 



4l| 



At first we introduce an effective flux as 

, 1-Sw 



(11) 



2 

As evident from Eq. ([3]) the spectrum is symmetric with respect to xs = 4>s/4>o- Further we 
define the shift As = xs — 1(1'), where S =f or L Here 1(1') = [x-j-m] where [x] means the 
nearest integer which is less than x. As is shown in Fig. 1. When A$ = 1, it is equivalent 
to A s = 0. Furthermore, A s = [1/2 — A s \ is the distance between the xs and the nearest 
half integer. 
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A. Spinless Particles 



For N spinless particles we distinguish two cases: N is an even or an odd integer. 
Case (1): If N is an even integer then 

4?en(A) = sgn(A)iV(i-A), 

4ven( A ) = ai« 2 ^(^r 1 )sin(A^6 r )cos(l - 2A)6 T , (12) 

where sgn(x) equals +1, for x > 0; for x = 0, and -1 for x < 0. Case (2): If N is an odd 
integer then 

J odd( A ) = ~ *) N (\ - A )> 

J« d (A) = a l a 2 h{Sl Tl )sm{Nb T )cos{l-2A)b T . (13) 

B. Particles with 1/2 spin 

For spin 1/2 particles we consider four cases. 

Case (0): For an even number of particles' pairs, i.e. N = 4m, where m is an integer we 
find 

4 0) (A 5 ) = 4?en(A s ), lg\A s ) = I^en(^s), (14) 

Case (1): For an odd number of particles' pairs, i.e. N = Am + 2 we obtain 

4 0) (A 5 ) = /^ d (A 5 ),7«(A 5 ) = /« d (A 5 ). (15) 

Case (2): For an even number of pairs plus one extra particle, i.e. N = 4m + 1 (there is 
one particle whose spin is unpaired as compared with case (0)) we find 

J ext,5( A ^) = -sgn(A s )sgn(i-A 5 )(^ + i-A 5 ), 

4xt,s( A s) = a ia2 h(n Tl ) sm(b T ) cos(^^ + 1 - 2A 5 )6 r . (16) 

To determine which spin state is occupied by the extra particle one compares the distance 
of the symmetric axis to the nearest half integral axis, i.e. A s . The one with the larger 
distance will be occupied. 

Case (3): For an odd number of pairs plus one extra particle, i.e. iV = 4m + 3. Here we 
use case (1) and determine the contribution to the current from the extra particle 

Cs^) = ^( A ^gn(i-A s )(^ + i + A 5 ), 

7 



Which spin state is occupied by the extra particle is governed by A s . The level with the 
smaller As is populated. 

V. SPIN-DEPENDENT CURRENT (SC) 

In presence of a static magnetic field and the SOI but in the absence of the pulse field 
the PCC is accompanied with a persistent SC (PSC). Switching on the pulse field generates 
a spin-dependent charge currents due to the SOI, and also a dynamic SC (DSC) that can 
be controlled by the parameters of the pulse field. The SC density is 

j:o,5o(r^^) = K{^S• t (r / ,^)v's^(r^^)} ) 

where 

s = (h/2)d z 5(r' - r) 
is the local spin density. The SC associated with level n , So i s 

1 /"27T PT2 

IS no,s«(t>T) = ^J o d V J^ dr% o ^{v\t) (18) 
and can be evaluated as 

=I^ElQ 5o (no,So,t)| 2 J Df«, (19) 

n 

where 

I s0 = e ip E a/(2n) 

sets the unit SC and 

Dt = i(a So ) 2 ~ (b s °mn --£-)- (^) 2 . (20) 

00 

Here 

(a 5 °) 2 -(6 5 °) 2 = S cos 7 , 
and So = ±1. The SC after applying two pulses to the ring is a sum of two parts 

Kl,s (t >r) = I£j£(f > r) + r-f o (t > r), (21) 



8 



where 



K:f (t>r) = I s o[^cos T ][(n-f) + 

<Po 



S 



I s0 [5 , cos7] 



AO) 

no, So 



2 cos 7 



(f>r) 



is the static PSC [27j and the DSC part is 



C,sb( t>r ) = ho[S cos^][a 2 (cos<f) no , So {r)], 



- I^cost]-^" 
Summing over all occupied energy levels we find 



I?}e.(t>T) 



l S s (t>r)=r s f\t>r)+r s f\t>r) 



S 



where (4° n ),(1) (t > r) are PCC and DCC) 



I S sf Ul \t>r) = I s0 [S cos7] 



j£>'«(t>r) 



(22) 



(23) 



(24) 



(25) 



VI. NUMERICAL RESULTS AND DISCUSSIONS 



We performed calculations for pulse-driven ballistic quantum rings fabricated by an ap- 
propriate confinement in a quantum well of In x .Ga!_ x As/InP [42]. Our results are also valid 
for other III-V or II- VI semiconductor quantum rings with spin orbit, e.g. GaAs-AlGaAs 

n 

quantum well, or HgTe/HgCdTe quantum ring [43J. We shall present the total charge current 
(TCC) which is a sum of PCC and DCC over all the occupied states. Total spin-dependent 
current (TSC) is obtained in the same way [281 ]. 

Fig. 2 shows how the flux and the SOI affect the PCC, DCC and TCC. Without the 
SOI, the jump of the PCC occurs at integer flux for even pair occupation, shown in Fig. 
2. The jumps are different in other occupations (see 28J), here we only focus on the even 
pair occupation case for clarity. The periodic sawtooth dependence of the PCC on the 
flux exhibits has been studied before, e.g. 27|, |28j. At finite SOI the jumps in PCC are 
shifted to 0/0o = I + (1 T w )/2; the two jumps are the consequences of a superposition 
of the contributions from the two spin channels. When the SOI strength is such that 
7 = — arccos(l/2n), (n = 1, 2, • • • ) the two jumps become at the half integer which is just 
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FIG. 2: (Color online) Persistent charge current (PCC), dynamic charge current (DCC) and total 
charge current (TCC) are shown in (a), (b) and (c) respectively. The spin orbit angle are 7 = 
0°,— 20°,— 40° and —50° for the solid lines, the dash lines, the dot lines and the dash-dot lines 
respectively in all graphs. The other parameters are N = 100, a = 100 nm, r = 26.3 ps, F\ = 
F 2 = l kV/cm. 

the case of An + 2 occupation in absence of the SOI [28]. The slope ratio between the two 
jumps is the same. As can be inferred from the analytic expressions DCC (cf. Fig. 2) 
depends smoothly on the flux. SOI results in a phase shift moving or even exchanging the 
positions of the minima and maxima, as is for 7 = —60°. The origin of the shape of TCC 
is deduced from those of PCC and DCC. Here the magnitudes of the two contributions is 
crucial: The PCC magnitude is related to the numbers of charge carriers, while the DCC 
magnitude is determined primarily by the product of the a± and a 2 (that can be externally 
varied by changing the pulse intensities), the delay time r, and the ring radius. 

Fig. 3 shows the TCC dependence on the ring radius (in the absence of the SOI). As 
expected, a larger a enhances the DCC. On the other hand, a enters the Bessel function 
argument whose increase suppresses the magnitude of DCC. It can be shown that the period 
of the oscillation with r increases with increasing the radius. The magnitudes of the maxima 
and minima are larger with larger radius. 



Now we discuss the spin-dependent current projected onto the z direction [271 ]. i.e. I Sz . 



The spin-dependent current projected onto the 7 direction (e.g. the quantization axis of the 
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FIG. 3: (Color online) Contour plots of TCC on (f> and r are shown for different radius of the 
ring, a = 100 nm, 200 nm, 300 nm and 400 nm in (a), (b), (c) and (d) respectively. The other 
parameters are N = 100, 7 = 0°, F\ = F 2 = 500 V/cm. 



local spin frame) is I Sz = J 7 cos 7 28]. 




FIG. 4: (Color online) PCC and PSC, DCC and DSC, and TCC and TSC are shown in (a), (b) 
and (c) respectively. The solid lines are the charge current, and the dash lines are the spin current. 
7 = —40°, the other parameters are the same to Fig. 2. 

PSC posses steps at PCC jumps (Fig. 4(a)) as a function of (p. This can be understood 
from the ratio of PCC for different spins; the SOI only introduces a relative effective flux 
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shift (see Eq. ( ITT!) ) leading to a constant spin-dependent current between the jumps. For 
DSC vs. cf) (Fig. 4(b)) the effective flux leads to a shift of the DSC along the flux axis. The 
physics behind this shift is that SOI provides a SU(2) flux, meaning that the pulse-driven 
(local frame) spin-up electrons experience a different flux than those with down spin, leading 
to a substantial spin-dependent current. In contrast, a static magnetic flux does not induce 
a spin-dependent current in the absence of the SOI. This shift and the jumps in the step 
function of the PSC explain the behaviour of TCC in Fig. 4(c). 




FIG. 5: (Color online) TCC (a) and TSC (b) vs. pulses strengths F\ = F2 and delay time r. 
7 = —40°, 4> = 0.2. Other parameters are ./V = 100, a = 400 nm. 

The control of TCC and TSC by tuning the pulse-field parameters is demonstrated in Fig. 
5. Because the two pulses transfer a net angular momentum setting the electrons in motion 
but they do not couple directly to the spins the TCC and TSC show the same pattern with 
r and F. From an experimental point of view it is essential to note that we are dealing with 
non-equilibrium quantities which opens the way for their detection via their emission. E.g., 
the TCC and TSC can be detected by measuring the current-induced magnetization of the 



ring and the generated electrostatic potential 



44 1 



VII. GENERAL AND CONCLUDING REMARKS 

In presence of spin-orbit coupling, two time-delayed appropriately shaped electromag- 
netic pulses generate spin-dependent charge currents. As shown previously for the spin- 
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independent case 171 ]. the sign of the current and its magnitude are controllable via the 
delay time and the strengths of the pulses. From a symmetry viewpoint, similar phenomena 
may be expected to occur for other geometries (wires, squares, etc). However, as shown for 
unbiased superlattices 45( (without SOI) details of the generated currents may differ quali- 
tatively. Application of an appropriate train of pulses open the possibility of controlling or 
even stopping the current [4l|. For increasing the magnitude of the current more intense 
pulses should be applied. 

For generating currents in quantum rings one may also apply circular polarized laser 
pulses fisl . \v\ . In this context we note the following: From an electrodynamics point 
of view, generating currents by our pulses is a completely classical effect, i.e. currents 
are generated in a completely classical system, even though in our case the subsequent 
excited carrier evolution is quantum mechanical. For this reason our current is robust 
to disorder and geometry modifications. In addition, the tunable time delay between the 
pulses allows an ultra-fast control the current properties. Using circular polarized laser pulses 
generates currents for quantized systems (in which case the rotating-wave approximation can 
be applied) . For systems with level broadening on the order of level spacing no appreciable 
current is generated. Our disadvantage however is that our pulses are much more demanding 
to realize experimentally, whereas laser pulses are readily available, in particular at high light 
intensities allowing thus for a strong current generation. 

The DSC is proportional to the DCC which can be comparable to the PCC for small or 
moderate occupation number case as seen in Fig. (4b). The DCC depends on the strength 
of the field and the delayed time between the two pulses sensitively and dramatically. We 
provide now an explicit calculation for the typical values of the CC and SC. For In^Gai^As 
/InP quantum well 42j we have m* = 0.037mo- For a ring with radius of 100 nm, the line 
velocity is then about 5000m/s, and the current unit is I^/a ~ 8nA which corresponds to 
the angular velocity current for one particle. If we convert it into the unit of an induced 
magnetization it is a radius-independent quantity Mo ~2 meV/T (here we use the formula 
Mq 7ra 2 (/o/«) valid for rings considered here 17|). 



The work is support by the cluster of excellence " Nanostructured Materials" of the state 
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